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1 Introduction 



In Ref. [T] (hereon referred to as (I)) a simple model, representative of QCD, 
was introduced and applied to the calculation of the spectrum of mesons. 
It is a Lipkin type model j2] for the quark sector, coupled to a boson level 
which is occupied by gluon pairs with spin zero. The four parameters of the 
model were adjusted in order to reproduce 13 known meson states with spin 
zero or one. The calculated spectra, for mesons with spin different from the 
ones used in the fit, were found to be in qualitative agreement with data. As 
reported in (I), the calculated meson states contain many quarks, antiquarks 
and gluons. The gluon contributions were found to be of the order of 30%. 
The model predictions (I) are free of the so-called multiplicity problem, i.e. 
that a given state can be described in many ways, which is removed due 
to the action of particle mixing interaction. The model itself resembles the 
one of Ref. jSj which treats nucleons coupled to pions. Also, it is related 
to the work of Ref. [3] , which describes quarks and uses particle conserving 
interactions. Generally speaking, the model of (I) belongs to the class of 
models described in Refs. El- The gluon part in (I) is fixed jTJ and 
does not contain any new parameters. The validity of the basic theoretical 
assumptions, and the applications to low and high temperature regimes, has 
been studied for mesons with flavor (0,0) and spin [8 . The aim of these 
studies was to formulate a manageable, schematic, albeit realistic, model 
to describe qualitatively QCD at low and high energies. Since the model 
is algebraic, i.e. all matrix elements are analytic, and exactly solvable, it 
can provide a non-perturbative description based on QCD relevant degrees 
of freedom, like quarks, antiquarks, and gluons. This, in turn, allows to 
test other microscopic many body techniques previously applied to the non- 
perturbative treatment of real QCD HQ] . Although the proposed model (I) 
is probably too simple to describe real QCD, it contains all basic ingredients 
of real QCD. These are the correct number of degrees of freedom associated to 
color, flavor and spin, and the orbital degree of freedom , which is contained 
in the degeneracy 2fl of each of the quark levels. 

In this work we investigate the behavior of the model, in the finite temper- 
ature regime. By starting from the model predictions of the meson spectrum, 
we calculate the partition function and different thermodynamical quantities, 
like the energy density and the heat capacity as a function of temperature. 
Next, we focus on the calculation of meson production rates. As we shall 
show, these production rates are in qualitative agreement with the experi- 
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ments. Also, we calculate absolute production rates for pions and kaons. Fi- 
nally, we concentrate on the transition from the Quark- Gluon-Plasma (QGP) 
[TH IT2"] to the hadron gas. The results support the notion that the present 
model may be taken as a toy model for QCD. 

The paper is organized as follows: In section 2 the model is shortly out- 
lined, since the details have been presented in (I). In section 3 we calculate 
the partition function and give the expressions for the relevant observables. 
In section 4 the model is applied to the description of the QGP. There, we 
present and discuss the results corresponding to some branching ratios and 
absolute production rates. Finally, conclusions are drawn in section 5. 

2 The model 

As described in (I), the fermion (quarks and antiquarks) sector of the model 
consists of two levels at energies +Uf and — ujf, each level with degeneracy 
2Q = n c rifn s , where n c = 3, rif = 3 and n s = 2 are the color, flavor and 
spin degrees of freedom, respectively (see Fig. 1 of (I)). Each level can be 
occupied by quarks. Antiquarks are described by holes in the lower level. 
Equivalently, one can use only the positive energy level and fill it with quarks 
and antiquarks with positive energy. The Dirac picture is useful because it 
gives the connection to the Lipkin model as used in nuclear physics. The 
quarks and antiquarks are coupled to gluon pairs with spin zero. The energy 
of the gluon level is 1.6 GeV [7j, and the energy u)f is fixed at the value 
ujf = 0.33 GeV , which is the effective mass of the constituent quarks. 

The basic dynamical constituent blocks of the model are quark-antiquark 
pairs B X f SM which are obtained by the coupling of a quark and an antiquark 
to flavor A (A = 0, 1) and spin S (S = 0, 1). The index / is a short hand 
notation for hypercharge Y, isospin I and its third component I z . Under 
complex conjugation the operator obeys the phase rule defined in [TB"] . 

The states of the Hilbert space can be classified according to the group 
chain 

[1"] [h] = [hh 2 h 3 ] [ft] 

U(4n) D 17(g) <g> 17(12) 

u u 

(A c ,/x c ) SU C (S) (A /)M/ ) SU f (3)®SU s (2)S,M , (1) 
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where the irreducible representation (irrep) of the different unitary groups 
are attached to the symbols of the groups. The irrep of U (40) is completely 
antisymmetric (fermions) and the one of U(^), the color group U(3) for 
0=9, and C/ (12) are complementary [TI] . The color irrep (Ac/ic) of the 
color group SUc(3) is related to the hk via Ac = hi — h 2 and \ic = h 2 — h 3 . 
The reduction of the U(12) group to the flavor (SUf(3)) and spin group 
(SUs(2)) is obtained by using the procedure described in [TO] . In (JTJ) no 
multiplicity labels are indicated (see (I)). 

The classification appearing in (JTJ is useful to determine the dimension 
and content of the Hilbert space. Instead of working in the fermion space 
we have introduced a boson mapping [T71 ITS] . The quark-antiquark boson 
operators are mapped to 



a Xf,SM ~^ u Xf,SM 

B\f,SM — ► b X f,sM ■ (2) 

where the operators on the right hand side satisfy exact boson commutation 
relations. 

The model Hamiltonian is defined completely in the boson space and it 
is given by 



// 



2ufUf + uj b n h + 



XS 



2b[ s b xs 
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)b+ 



[b[ s r+2b[ s b xs+ (b xs ) 



(3) 



where (b\ s ) 2 = (b\ s ■ b\ s ) is a short hand notation for the scalar product. 

Similarly for {b\s) 2 and (b\ s b xs ). The factors (1 — ) simulate the terms 
which would appear in the exact boson mapping of the quark-antiquark pairs. 
The and b are boson creation and annihilation operators of the gluon pairs 
with spin S = and color A = 0. The interaction describes scattering and 
vacuum fluctuation terms of fermion and gluon pairs. The strength V\g is 
the same for each allowed value of A and S, due to symmetry reasons, as 
shown in (I). The matrix elements of the Hamilton operator are calculated 
in a seniority basis. The interaction does not contain terms which distinguish 
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between states of different hypercharge and isospin. It does not contain flavor 
mixing terms, either. The procedure used to adjust the four parameters 
(values of V\s), was discussed in detail in (I). 

The disadvantage posed by working in the boson space is the appearance 
of un-physical states. In (I) we have presented a method which is very efficient 
to eliminate spurious states, as we shall show in this paper. The suitability 
of the Hamiltonian (jSJ) to describe the gluon pair and quark- ant iquark pair 
contents of mesons has been discussed in details in (I). 

As a next step, in this paper, we have introduced temperature and dis- 
cussed the transition to and from the QGP. As it can be expected, because 
of the schematic nature of the model, we may attempt to describe only the 
general trends of the observables. To achieve this goal, further assumptions 
have to be made respect to the volume of the system because the model, as 
it has been proposed in (I), has no a priory information about the volume of 
the particle. 

3 The partition function, some state variables 
and observables 

The group classification of the basis ((TJ), allows for a complete book-keeping of 
all possible states belonging to the Hilbert space of (jBJ). The corresponding 
partition function, which contains the contribution of the quark- ant iquark 
and gluon pairs configurations introduced in the previous section, is given by 

Z qgl = $>m(Ac,A*c) E £(2.7+1) ££ 

[h] (A /( a*/) J P=± i 

multiE^e-^-^ 13 -^-^ , (4) 

where [is, Us and are the baryon, strange and isospin chemical potentials, 
respectively. The sum over [h] = \h\h2h3\ denotes all color irreps of U(3) with 
J2k hk = N, where N is the total number of quarks in the two levels (Dirac's 
picture). The transposed Young diagram [hf, obtained by interchanging rows 
and columns, denotes the U(12) irrep. The index "i" refers to all states with 
the same color, flavor, spin and parity (P). These states are obtained after 
the diagonalization of the Hamiltonian (J3J). For mesons belonging to the tx-t] 
and io-p octet, the mass values entering in (jlj) do not take into account flavor 
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mixing. The eigenvalues £7, are denoted by the eigenvalue index i, and they 
are functions of all the numbers needed to specify the allowed configurations, 
namely: s, (A/, fif), P, S, [h] and of the cutoff for the different boson species 
[A, S] (see (I)). The quantities B and I z are the baryon number and the third 
component of the isospin. According to the experimental evidences the value 
III — is a reasonable approximation, and we have consistently adopted it 
in our calculations. The dimension corresponding to color configurations is 
given by mult(Xf,fif) = |(A/ + !)(/// + 1)(A/ + /// + 2) and of the spin by 



Since the eigenstates of the Hamiltonian Q have been calculated after 
performing a boson mapping, as described in (I), we have consistently fixed 
the corresponding cutoff- values at the values 2Q, Q, 4p, and ^ for the boson 
pair species [0,0], [0,1], [1,0] and [1,1] respectively (see (I)). These values are 
adequate when the fermion (quark-antiquark) configurations entering in the 
boson states correspond to a full occupation of the fermion lower state (—Uf). 
These numbers may be modified when the fermion configurations correspond 
to states where the upper level is partially occupied and the lower level is 
partially unoccupied. The distribution of the occupation in the upper and 
lower fermion levels is fixed for a lowest weight state \lw > of a given U(12) 
irrep, defined by B\f tS M\lw >= 0. The irrep of U(12) is given by a Young 
diagram ^3] with rrik boxes in the fc'th row. The lowest weight state is given 
by J2t=i m k quarks in the lower level and Y^k=i rn k quarks in the upper level. 
The highest weight state is obtained by interchanging the occupation. The 
difference of the number of quarks in the upper level, appearing in the highest 
and lowest weight states, gives the maximal number of quarks we can excite 
for a given U(12) irrep. This number is given by 



For the case [3 6 6 ], used in (I), we found J — Q. Therefore, 2 J is the 
maximal number of quarks we can shift to the higher level, i.e. it is equal 
to the maximal number of quark-antiquark pairs which can be put on top of 
the lowest weight state of a given U(12) irrep, which is also the state with 
the lowest energy in absence of interactions. 
The total partition function is given by 



(2J+ 1). 



6 12 



2 J = m k - m k 



(5) 



k=l k=7 



Z 




(6) 
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where Z g is the contribution of all gluon states [7] which do not include con- 
tributions of gluon pairs with spin zero. It is written Z g = J2a ex P(~PE a ). 
The values E a can be deduced by using Eq. (40) of Ref. [7j. Except for 
the gluon pairs with spin zero, all other gluon states are treated as specta- 
tors because the Hamiltonian in © does not contain interactions with these 
other states. Note, that the interaction between the gluons is taken into 
account explicitly in the model of Ref. [7j. As a short hand notation we will 
abbreviate the partition function by Z = J2i e - ^, taking into account that 
Si contains the information about the chemical potential, the contributions 
of the quarks and gluons. 

The observables < O > are calculated via [TT?] 

<0>= ^ Ue , (7) 

where the index a denotes the color configurations, i.e; a=(0,0) when only 
color zero states are considered and a = c when also states with definite color 
are allowed. This distinction is needed to investigate the phase where color 
confinement is effective and the phase where color is allowed over a wide area 
of space. The quantities to calculate are the internal energy (< E >), heat 
capacity (< C >), average baryon number (< B >), strangeness (< s >) and 
the expectation value of different particle species (< n& >), where k refers to 
the quantum numbers of a particular particle and J2k < n k >= 1, with the 
sum over all possible quantum numbers. 

The particle expectation values have a simple expression because they 
select one of the eigenvalues at the time, thus if the state of a given particle 
is denoted by "z" and E^ is its energy, the particle expectation value is given 
by 



e -f3(E t -fi B B-fM 3 s) 

<rii> = , (8) 

(where we have used the value /i/ = 0, for the isospin chemical potential ). 

At this point we have to make an assumption upon the volume consid- 
ered. The whole reaction volume can be divided in elementary volumes, and 
we assume that the elementary volume (V e i) is of the size of a hadron, cor- 
responding to a sphere with a radius of the order of 1 fm. Later on we shall 
show that this choice is reasonable, as seen from the calculated thermody- 
namic properties of the whole system. Another assumption is related to the 
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interaction, which does not take into account confinement. We shall discuss 
two scenarios, namely: a) no additional interaction related to color is taken 
into account for temperatures above a critical (de-confinement) value, and b) 
confinement is operative for temperatures below the critical de- confinement 
temperature. In the regime (a) the lowest state with color (1,0) lies at the 
energy 2a;/, and it corresponds to put one quark in the upper fermion level. 
Although it is a possible configuration, the Hamiltonian (j3J) cannot act upon 
it. This is consistent with the fact that above a certain temperature, T c 
(de-confinement temperature), only color non-singlet states are allowed. The 
actual value of T c will then give us an idea about the regime where hadroniza- 
tion is operative. In the real world hadronization, i.e. confinement, should 
set in below a critical temperature, as a true phase transition. In our model 
this will be signaled by a sharp transition from a state where color non-singlet 
states are still allowed (T > T c ) and a state where confinement is effective 
(T < T c ). 

4 Description of the high temperature regime: 
the QGP 

We shall first discuss the case where no additional color interaction is 
taken into account. The states with energy Ei for a given flavor, spin and 
parity are obtained from the diagonalization of the model Hamiltonian, now 
including flavor mixing and the corrections due to the Gel'man-Okubo mass 
formula for the two lowest meson nonets (one with spin zero and the other 
with spin 1). In Figure 1 we show the internal energy as a function of 
the temperature T, with and without interactions. The results shown by 
a dashed line have been obtained by calculating the internal energy in the 
fermion space without interactions. In this case the Hamiltonian has a simple 
image in the fermion space and the calculation of the partition function can 
be performed exactly. The dotted line shows the results obtained by working 
with the boson mapping, and by enforcing the corresponding cut-offs in the 
maximal number of bosons as explained above and in (I). The internal energy 
is a good indicator of the number of active states in the Hilbert space. Note 
that the results shown by both curves, the doted and dashed lines, practically 
coincide. This suggests that the number of active states is nearly the same 
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in the boson and fermion spaces for a wide range of temperatures. This does 
not imply that all un-physical states have disappeared but rather that the 
approximate method of cutting un-physical states works reasonable well. The 
curve shown by a solid line, in the same Figure 1, gives the internal energy 
obtained from the calculation performed in the boson space and in presence 
of interactions. Although the curve does not show a clear phase transition 
of first order (e.g: a sharp increase of the energy in a narrow interval around 
T c ) the behavior around T c = 0.170 GeV is pretty suggestive of it. We have 
interpreted the observed smearing-out of the curve as follows: for T = 
the vacuum state is dominated by pairs of the type [1,0] and in this channel 
a quantum phase transition [S] does indeed take place. By this we mean 
that the pairs [1,0] are effectively blocked at high temperature. The other 
channels contribute less significantly to the ground state (see (I)) and interact 
weakly than the ([1,0]) channels, therefore, they remain in a perturbative 
regime. Thus, as the temperature increases, an approximate first order phase 
transition takes place in the channel [1,0], a mechanism similar to the one 
shown in Ref. [S] for the [0,0] channel, while the other configurations remain 
un-affected. The superposition of these two mechanism leads to the smearing 
out of the curve around T c , as shown in Figure 1. 

In Figure 2 the expectation value of the Casimir operator (C 2 =< C2 >) 

of color and its variation (AC 2 = \J < C 2 2 > — < C 2 > 2 ) are shown. The 
eigenvalue of the Casimir operator, for an irrep with color numbers (Ac Hc)i 
is given by C 2 (Xc, ^c) — + ^c^c + ^c* + 3(Ac + He)- As a reference, for 
a color (1,0) irrep C 2 = 3 while the irrep (1,1) has C 2 = 9. We assume that 
fiB and n s are zero. It is interesting to observe that, in the present model, the 
variation of the color is approximately symmetric around T = 0.170 GeV. 
A possible interpretation is the following: at high energy the probability 
to have a color non-singlet state is large (the variation is not large enough 
to allow color singlet states) and a QGP is formed where color is effective 
over a wide range in space. From T = 0.170 GeV on the probability to 
find a state in color (0,0) is significantly increased, since the variation is 
large enough to allow color singlet states. In lowering the temperature the 
variation is much larger than the average color and the whole QGP dissolves 
in droplets of color zero. Within the present model, these results, of the 
average color and its variation, are signals of the transition to the hadronic 
phase. Accordingly, we assume that it takes place near T c = 0.170 GeV, 
for jiB = Us = 0. We may now calculate the bag pressure and construct 
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Figure 1: Internal energy as a function of T. The dashed line corresponds 
to the exact result in the fermion space, without interaction. The dotted line 
is the internal energy calculated in the boson space with the cut-off for the 
different boson species, as discussed in the text. The solid line represents 
the results obtained with the full Hamiltonian and calculated in the boson 
space. No additional color interaction is taken into account. 
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Figure 2: The expectation value of color (C2 =< C2 >, solid line) and 

its variation (AC2 = \j< C 2 2 > - < C 2 > 2 , dotted line) as a function of 
T. The variation reaches the same value as the expectation value around 
T = 0.170 GeV. The values /j,b = Us = were used in the calculations. 
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the T — fi B diagram. At T = 0.170 GeV the pressure is determined via the 
expression p = Tl ™( z ) ; where $ = —Tln{Z) is the grand canonical partition 
function [TH] and V e i is the elementary volume V e i = -f-r e z- For r e i=l fm we 
obtain a bag pressure p^ of about 0.17 GeV, which is in reasonable agreement 
with standard values. When the chemical potentials \x B and p, s are different 
from zero, the temperature dependence of the internal energy changes and 
also changes the value of the temperature for which the pressure is equal to 
the bag pressure p± = 0.18 GeV. The results are shown in Figure 3. Assuming 
that the local strangeness is < s >= 0, we arrived at a functional relation 
between T, /i B and /i s , i.e. f(T, fi B , /-0 = 0, which fixes function of 

fi B - The results of this functional relation are displayed in Figure 4. 

Once the chemical potential /x s is adjusted, by using the results shown in 
Figures 3 and 4, the chemical potential /i B and the transition temperature 
T c can be consistently determined. 

Up to now we did not take into account an interaction which generates 
confinement. This has to be done by hand. One possibility is to assume 
that the transition from the QGP to the hadronic phase takes place within 
a very small range of temperatures around the critical temperature T c . We 
require that the partition function above T c allows any color while for T <T C 
it contains only color zero states. Finally, chemical equilibrium connecting 
both phases, the QGP and the hadron gas, is understood. Figure 5 shows the 
results of the internal energy, without confinement (upper curve) and with 
confinement (lower curve). The solid line connecting both curves indicates 
the values for which confinement vanishes above T c . As seen from the results, 
the transition is now of first order. Figure 6 shows the heat capacity calcu- 
lated for the case without confinement (dotted line) and with confinement 
below T c (dashed line). The solid line interpolates between them, as in the 
case of the internal energy (see Figure 5). 

The model can first be tested in the energy region below the transition 
temperature T c , where the hadron gas should prevail. The confinement is 
effective and therefore we have to use the partition function ^ a= (o,o) i n the 
equations (J7J) and (jHJ). We take, as an example, the measured total production 
rates of ir + at 10 GeV/A, as reported in the SIS-GSI experiment The 
system considered was Au+Au and we assume that all particle participate, 
i.e. iVp or t=394. In Figure 2.3 of Ref. [20] a temperature of about T = 0.13 
GeV is reported. Assuming local strangeness conservation , < s >= 0, we 
obtain a relation of fi s versus T, which is depicted in Figure 7. From there 
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Figure 3: Temperature versus chemical potential. The curves gives the 
boundary where the pressure of the system in an elementary volume of radius 
r = 1 fm is equal to the bag pressure. The values of the temperature T are 
given as a function of the chemical potentials /j,b (solid line) and fi s (dashed 
line). The label \x denotes both /i B and ji s . 
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H (GeV) 



Figure 4: Chemical potential n s as a function of the chemical potential //£, 
for zero local strangeness < s >=0. To each pair of values {hb and fi s ) it 
corresponds a temperature T, which is determined from the limiting values 
shown in Figure 3. 
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Figure 5: Internal energy. The upper curve shows the results obtained 
without considering confinement. The lower curve corresponds to the results 
obtained with the inclusion of confinement. The solid line shows the results 
obtained with the partition function where confinement is considered to be 
fully operative below T c (//b=// s =0). 
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Figure 6: Heat capacity of the system. The solid line interpolates be- 
tween the case with confinement below T c to the case without it above T c 
(a*b=A*s=0). The dotted and dashed line correspond respectively to the cases 
where color is allowed and without color. 
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we obtain, for the reported temperature, fi s ~ 0.128 GeV and via Figure 4 a 
value hb ~ 0.55 GeV. In Figure 8 we show for a fixed value of physically 
acceptable at temperatures near T = 0.13 GeV, the resulting total production 
rate of tt + (N n ) as a function in the temperature T. For T = 0.13 GeV the 
production rate is approximately 180 pions tt + , which is close to the value 
160, which we have obtained by using Figure 2.3 of Ref. [20j. The good 
qualitative agreement with the experiment demonstrates that the present 
model is able, indeed, to describe, approximately, observed QCD features. 

We have also determined ratios of particle production and some absolute 
production rates. The particle production is calculated for temperatures 
just below T c , where only color zero states are allowed. This implies that 
the partition function to use is Z(o,o)- We then apply Eq. (jSJ). Note that in 
the expression of the particle-production ratios the partition function cancels 
out and only the dependence on the mass of the particles and the chemical 
potential remains. Figure 9 shows results for some particle-production ratios 
for beam energies y/s = 130 A GeV. The experimental values are taken from 
Ref. [21], based on the experiment described in Ref. [22] (see also j2H])- For 
baryons, only the ratios of particle and anti-particle production are shown 
because these expressions are independent of the mass of the baryon. As 
noted in (I) the masses of the baryons are not well reproduced because they 
are considered as consisting of three idealized fermions on top of the meson 
sea. The interaction to the meson sea is not taken into account yet, but 
indications about how to do it are given in (I). 

The central value of the K~ / K + production ratio, shown in Figure 9, 
was reproduced with the values /i s = 0.012 GeV and \i# = 0.044 GeV. The 
other ratios are predicted by the model. Considering the simplicity of the 
model, the ratios are found to be in a reasonable agreement with data. 

In order to obtain the total yields for kaons and for the 7r + pions it is 
necessary to introduce further assumptions about the size of the QGP. The 
baryon density is given by ^7^, where V e i is the size of the representative 
volume, as explained before. In order to conserve, on the average, the baryon 
number we multiply the baryon density by the total volume and require 
that it must be equal to the total baryon number, given by the number of 
participants N part . This leads to the total volume 

where the index a refers to, as in the partition function, color (a = c) when the 
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Figure 7: The dependence of /i s on the temperature T, for /!# = 0.55 GeV 
and assuming < s >=0, for the SIS-GSI experiment (2Hj, Au+Au at 10 
GeV/A. 
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Figure 8: Total production rate of ir + for Au+Au at 10 GeV/A [201 - F° r 
T = 0.13 GeV the total n + production rate is approximately 180. The value 
fiB = 0.55 GeV was used in the calculations, as discussed in the text. 
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Figure 9: Some particle production ratios for the beam energies ^/s = 130 A 
GeV taken from Ref. [21]. The calculated values are shown with full squares 
and the experimental values are shown within error bars. 
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average value is calculated in the QGP and a = (0, 0) when it is calculated in 
the hadron gas. In the QGP the average value < B > c before the transition 
is smaller than the average value < B >(o,o) m the hadron gas after the 
transition. This is due to the small value of Z(o,o) a t T c , since many other 
possible color states are excluded from it which do contribute to Z c . As a 
consequence, for T = T c the volume of the QGP phase, as a function of N part , 
is much smaller than the one in the hadron gas phase. Assuming a sphere, 
the radius of the QGP phase is about 8 fm, and it changes to about 20 fm 
after the transition. This implies volumes of the order of approximately 2 
10 3 fm 3 and 3.4 10 4 fm 3 , respectively. 

This transition is, as pointed out earlier, assumed to take place suddenly 
at T c (probably it should be smeared out but we cannot describe it with 
the present model). This implies that within the scenario assumed there is 
a rapid expansion of the volume caused by the transition from the QGP to 
the hadron phase, which should be observed as a large outward motion. One 
possible interpretation is that most of the pions are produced during the 
transition, liberating energy and provoking a rapid expansion of the system. 
The energy gained is represented by the jump between the lower and upper 
curves of Figure 5, but its origin cannot be explained by the present model, 
where confinement was shifted by hand. 

Figure 10 shows the total pion yield as a function of the temperature T, 
corresponding to the Au+Au collision. The upper curve describes the total 
yield when all nucleons participate, while for the lower one we have taken 
Npart = 250. This value agrees better with the experiment, as seen from 
the results, and it means that in the collision about 250 nucleons partici- 
pate in the QGP. In Figure 11 the total kaon production rate is displayed. 
The upper curve corresponds to K + and the lower one to K~ absolute pro- 
duction rates, respectively. In both cases N part = 250 was used, the same 
value used previously in the calculation of the pion yield. The ratio of the 
curves was already adjusted at the point corresponding to the K + / K~ ra- 
tio. The absolute production rate and the shape of the curve, however, is 
a prediction of the model (as far as we can talk about "predictions" within 
this toy model). Considering the simplicity of the model it is surprising that 
the absolute production rate is well reproduced. This feature is common to 
other thermodynamical descriptions of the transition from the QGP to the 
hadron-gas |2lH |2*T]. 

Finally, in Figure 12, we show the calculated expectation values of the 
number of quark and gluon pairs as a function of the temperature T. At T = 
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Figure 10: Total production rate of ir + . The upper curve is for N part = 394 
(Au+Au) and the lower one is for N part = 250. Data are taken from [21]. 
Because of \ii = 0, the production rate for n~ coincides with the production 
rate for 7r + . 
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Figure 11: Total production rate of K + (upper curve) and K (lower curve) 
for Np ar t = 250. Data are taken from [21]. 
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GeV the results correspond to the fractions of gluon pairs and fermion (quark- 
antiquark) pairs in the physical vacuum state. At high temperatures the 
gluon part increases and takes over the fermion part, which shows saturation. 
However, at the temperatures of interest, i.e. around the point of the phase 
transition T c « 0.16 GeV, the gluon number is still suppressed with respect 
to the fermion pair number. This might be in favor of the ALCOR model 
which supposes a suppression of gluons in the QGP and takes only constituent 
quarks and antiquarks into account. Note, that at T = 0.170 GeV still a 
sensible amount of gluon pairs are present. 

5 Conclusions 

We have presented a toy model of QCD. The model is described in (I) and 
in this paper we have focused on the thermodynamic properties, at equilib- 
rium, emerging from the model. We have calculated the partition function 
with and without color, and studied the temperature dependence of some 
observables, like the internal energy, the heat capacity, and the production 
rates of particles. The parameters of the model were determined in (I), ad- 
justing the meson spectrum at low energy. Without further parameters the 
internal energy, the heat capacity and some particle ratios were determined, 
as explained in the text. 

We have applied the model to the case of the Au+Au collision at 10 
GeV/A J2U] and shown that it can reproduce qualitatively the absolute pro- 
duction rate of 7r + . At this energy the QGP has not yet formed and, therefore, 
the results show that the model can be applied to study schematically the 
thermodynamics of a hadron gas. 

Next, we have applied the model to energies where one assumes that the 
QGP has been already formed. The absolute production rate of tt + and 
kaons were calculated, just below the transition temperature, by taking the 
number of participant nucleons (N part ) as an input. The agreement between 
calculated and experimental values was found to be satisfactory. Also, the 
resulting production rate was described reasonable well, once the /j, s chemi- 
cal potential was fixed to yield the correct (observed central value) ratio. 
Some mass-independent baryon-antibaryon ratios, were qualitatively repro- 
duced by the model predictions. 

This demonstrates that the model is able to describe the general trend of 
QCD, in the finite temperature domain, and the transition to and from the 
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Figure 12: Expectation value of the fermion (quark-antiquark) and gluon 
pairs. The symbols Z tot and Z 00 , on each curve, indicate if color was al- 
lowed. The upper indices / and g on the Z refer to fermion and gluon pairs 
respectively. 



25 



quark gluon plasma. 
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